Electrons in 2D periodic magnetic and electric fields: Avoided band crossings, 
metal-insulator transitions, and changes in the spreading direction 



Bjorn Naundorf 

Max- Planck- Institut fur Stromungsforschung and Institut fur Nichtlineare Dynamik der Universitat Gottingen, 

Bunsenstr. 10, 37073 Gottingen, Germany 

Roland Ketzmerick 

Institut fur Theoretische Physik, Technische Universitat Dresden, 01062 Dresden, Germany and 
Max- Planck- Institut fur Stromungsforschung and Institut fur Nichtlineare Dynamik der Universitat Gottingen, 

Bunsenstr. 10, 37073 Gottingen, Germany 
(Dated: February 2, 2008) 

We study the energy spectrum and wave packet dynamics of electrons in a two-dimensional 
periodic potential subject to a perpendicular spatially modulated magnetic field. We observe avoided 
band crossings, metal-insulator transitions, and changes in the spreading direction of wave packets 
in two new situations: (i) for a purely magnetic modulation and (ii) in the limit of weak electric 
and magnetic modulations, where the Landau band coupling can be neglected. In the second case 
the corresponding classical system is no longer chaotic but integrable, whereby these phenomena 
are not expected. We argue that they are due to tunneling. 



I. INTRODUCTION 

Since the early days of quantum mechanics the problem 
of electrons in a periodic potential subject to a magnetic 
field has been studied Q. While the energy spectrum of 
an electron in a magnetic field is discrete, the spectrum 
of an electron in a periodic potential has continuous en- 
ergy bands. The combined system leads to fascinating 
commensurability effects, depending on the number of 
magnetic flux quanta per unit cell of the potential. These 
effects are visualized in the Hofstadter butterfly Q and 
are subject to intensive research. While many studies fo- 
cus on the case of uncoupled Landau bands, which can 
be described by the Harper equation jj, the full sys- 
tem exhibits many more features 0, 13- 13- S H El • For 
example, unusual sequences of the Hall conductance ap- 
pear |]Jj], which recently have been observed experimen- 
tally with lateral superlattices on semiconductor hetero- 
junctions |3~TI | . Another important feature of the Lan- 
dau band coupling is avoided band crossings, which oc- 
cur due to a classically chaotic limit and lead to metal- 
insulator transitions In systems with direction de- 
pendent modulation strengths they induce a counterin- 
tuitive change of the spreading direction of wave packets 
into the direction of strong modulation 0] . 

In the last few years attention has been paid to the case 
of 2D periodically modulated magnetic fields both in the- 

ory d ei m iijiiii m mm mnana 

and in experiment 28, 29, 3(J,|^l|. Consequently, the fol- 
lowing question arises: Do avoided band crossings, metal- 
insulator transitions, and the unusual ch ang e in transport 
direction found for electric modulation carry over to 
magnetically modulated systems? 

In this paper we derive a matrix equation for the energy 
spectrum of electrons in 2D periodic electric and mag- 
netic fields. We observe avoided band crossings, metal- 
insulator transitions, and changes in the spreading direc- 



tion of wave packets in two new situations: Firstly, we 
find these phenomena in the case of a purely magnetic 
modulation. Secondly, we study the limit of weak electric 
and magnetic modulations, where the Landau band cou- 
pling can be neglected. For the energy spectrum we de- 
rive a ID discrete Schrodinger equation with next-nearest 
neighbor coupling. Its classical counterpart is an inte- 
grable ID Hamiltonian, where avoided band crossings are 
not expected. We show, however, that avoided crossings 
exist and argue that they are due to tunneling. These 
avoided band crossings lead to metal-insulator transitions 
and changes in the spreading direction of wave packets 
in an unexpected regime. 

The paper is organized as follows: We proceed in 
Sec. II with the introduction of the system and show in 
Sec. Ill its spectrum and dynamics. In Sec. IV we study 
the limit of weak Landau band coupling. The correspond- 
ing classical limit is shown to be integrable in Sec. V 
where we also discuss the origin of the avoided band 
crossings. We finally summarize our results in Sec. VI. 



II. SYSTEM 

The dynamics of an electron with charge — e and mass 
to in a potential V(r) and a magnetic field B(r) is de- 
scribed by the Hamiltonian 



H = 



1 

2 TO 



p + eA(r) +V(r), 



(1) 



where the vector potential A(r) is given by B(r) = 
V x A(r). We are interested in a 2D potential and a 
perpendicular magnetic field, both periodic in the x and 
y directions with period a and b, respectively. Both the 
potential and the magnetic field can be decomposed into 
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Fourier series, 



We haven chosen the vector potential A(r) as, 



V(r) - V rtS e 2lTl{rx/a+sy/b \ (2) 

(r,s)^(0,0) 



B(r)=(B + B r ^e 2 * i{rx/a+sy/b) ] e z . (3) 

(r,s)^(0,0) 



A(r) = {0,B x,0) 



- E 



Br. 



27ri ^— ' (r/a) 2 



(s/b) 



_ e 2ni(rx/a+ S y/b) (- S /b, r / a,0) . 



(4) 



r 



such that it commutes with the momentum operator p, 
which reduces the following algebra. For vanishing mag- 
netic and electric modulation, B TtS = V r>s = 0, the eigen- 
functions of H are the well known Landau level eigen- 
functions, 



We note that while the set of rational magnetic fluxes 
4>/$o is of measure zero within the real numbers, it is 
dense within the real numbers and thus may approxi- 
mate any irrational $/$o- The eigenequations of the 
magneto-translation operators read 



{x,y\u, k) 



iy/b 



Kl 2 /b 



I 



(5) 



M x+pa \ip(k x ,ky)) 



M y + b \lP(k X ,ky)) = \^(k X ,ky)), 



(8) 
(9) 



with v being the index of the harmonic oscillator eigen- 
state <p v (z) — exp(— z 2 /2)H v (z) [^3, the phase K € K, 
and the magnetic length I 2 = hj (eB ). In the presence of 
a modulation, the states \u,n) are no longer eigenstates 
of H. As they constitute a complete orthonormal set, we 
use them as a convenient basis. 

The magnetic flux $ per unit cell of size ab in 
units of the magnetic flux quantum, <£>o = h/e, is 
given by §/&o = abBo/(h/e). If this ratio is ratio- 
nal, <£/$o = q/p, the Hamiltonian commutes with the 
magneto-translation operators 



M x+pa = e'» /|2 e rt 
M y+b = e be> y 



(6) 
(7) 



and there are common eigenstates of H, M x+pa , and 
M y+b . This will allow to treat the problem numerically. 



with the magnetic Bloch phases k x ,k y € [0,1). Their 
eigenvectors \ip(k x ,k y )) can be constructed from the 
chosen basis using M x -\- pa \v, k) = \v, n + 2nq) and 
M y+b \v,k) = e lK \v,k), yielding 
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\i>(k x ,k y )) = 22 2J a Ut n(k x ,k y )\v,2'K(n + ky)) 



with 



1/—0 n— — 00 



a = e~ 27rlk *a 



(10) 



(ii) 



Inserting \tp(k x , k y )) into the Schrodinger equation, 
H\ip(k x ,k y )) = E\ip(k x ,k y )), and projecting onto the 
basis states (v,2ir(n + k v )\ yields a set of equations for 
the coefficients a v<n) 



J 



1 



OO OO 



v + I J a^ n (k Xi k y ) + } j 

u'—0 S— — 00 



E 

&V' ,71 — S {kx •) ky ) = — QjV,n{kxiky}i 

RUf- 



(12) 



with the cyclotron frequency uj c — eB /m. This set can be reduced to values n = 0, 1, . . . , q— 1 by employing Eq. ltTT|l . 
which allows to write a v ' )n - s — e~ 27Tlk * L(™~ S )/9J a v i ,( n - s ) modgj with [^J denoting the largest integer that is smaller 
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than or equal to x. The coefficients f are given by 

1 q ^ ^ B r — , . - ■ 



rril/ ,V 



E 

r 

(r,a)/(0,0) 



huj r "'>" 2tt P 4i {r-r')r'a 2 + {s-s')s'/a 2 ~' " 

(r',s')/(0,0) 
(r,s)^(r',s') 



B r ,s/B 



2n(r 2 + s 2 /a 4 ) \ p 



v'D r J?_ 1 „(sa + ir/a) + \J v' + lD r ^ +1 v (—sa + ir/a) 



(13) 



Here we have introduced a = \fbfa and 



fc=0 



fc! (i/'-fc)!(i/-fc)!' 



where p = ^J^(ira — s/a). The same result has been 

derived for a — 1 in a different way in Ref. [T3 |. 

In the following, we restrict ourselves to the lowest 
Fourier components in the electric potential and the mag- 
netic field as well as identical spatial periods a = b such 
that, 



V(x, y) = hw c ( v x cos 



2ttx 27ry\ 
i h v v cos 



a J 



B(x,y) = B Q 1 + 6^ cos 



2irx 2irx\ 
i 1- o„ cos I 



a J 



(15) 



(16) 



i.e. in the notation of Eqs. 10 and ©, ZJ-i.o = £?i,o = 
and Vb,-i = Vb,i == fuo c v y /2. 



III. SPECTRUM AND DYNAMICS 

Figure H shows the energy spectrum for a purely mag- 
netic modulation, v x — v y — 0, and for different mag- 
netic modulation strengths in the x and y directions, 
b x /by — 0.8. One can see how the isolated Landau lev- 



els at b x 



are split up into a complex subband 



structure for b x + b y > 0. Upon further increase of the 
magnetic modulation strength, these bands broaden but 
typically repel and do not cross each other. This phe- 
nomenon has also been termed avoided band crossings 
[u3 |. In the left panel a subspectrum is shown where only 
k x is varied and k y = is kept fixed, whereas in the right 
part k y is varied, keeping k x = fixed. A clear duality 
of these subspectra can be seen with wide bands in one 
panel corresponding to narrow bands in the other panel. 
In the limiting case of irrational $/$o this corresponds 
to absolutely continuous and pure point spectral regions, 
respectively. Therefore, even in the rational case, where 
we present the numerics, we will denote the narrow bands 
as levels. 



Such spectral properties go along with directional 
properties of eigenstates, which was demonstrated for 
a purely electrostatic modulation in Ref. For ex- 

ample, a band under variation of k x and a level under 
variation of k y correspond to eigenfunctions extended in 
the x direction and localized in the y direction. This has 
immediate consequences for the dynamics of wave pack- 
ets. An initially localized wave packet, constructed from 
states in an energy interval with bands under variation of 
k x , will spread ballistically in the x direction and will lo- 
calize in the y direction. This is demonstrated in Fig. [2t- 
The initial wave packet is constructed from a symmetric 
eigenstate of the v — 2 Landau level that is projected 
on the subspace spanned by the eigenstates of the chosen 
energy range. 

In Fig. we present the spreading behavior for 
stronger magnetic modulation at the same ratio b x /b y — 
0.8 and for a similar energy range. We find that the 
spreading direction is switched to the opposite direc- 
tion. This nonintuitive dynamical effect was previously 
observed in the case of a purely electrostatic modula- 
tion It can be understood from the subspectra 
shown in Fig.^ Upon the increase of the magnetic mod- 
ulation strength b x +b y one observes metal-insulator tran- 
sitions in one subspectrum that go along with insulator- 
metal transitions in the dual subspectrum. In Ref. ^3] it 
was argued that such metal-insulator transitions can be 
induced by avoided band crossings, which are ubiquitous 
in the spectrum of Fig. 

Avoided band crossings and avoided level crossings are 
a generic phenomenon in quantum systems with a classi- 
cally chaotic limit. In fact, the present system with mag- 
netic modulation was shown to be chaotic [33] . Therefore 
the observed change in the transport direction by increas- 
ing the magnetic modulation strength can be understood 
along the line of reasoning previously applied to the case 
of a purely electrostatic modulation: A classically chaotic 
limit leads to avoided band crossings in extended sys- 
tems. They may induce metal-insulator transitions [l2| 
which lead to a change in transport direction . In ad- 
dition to these findings for a purely magnetic modulation 
we have verified that such changes in the transport di- 
rection just as well occur in the general case of combined 
electric and magnetic modulation. 

In the second part of the paper we want to demonstrate 
that there is a second origin of avoided band crossings 
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£>r + by b T + by 



FIG. 1: Spectrum of Eq. JjQl vs. the strength of the magnetic modulation, b x + b y , for a fixed ratio of modulation strengths, 
bx/by — 0.8, and without electrostatic modulation, v x = v y = 0. The average magnetic flux per unit cell is given by $/$o = 
55/89, which is an approximant of the golden mean (v5 — l)/2. In the left part, k x is varied, keeping k y — fixed, while in the 
right part k y is varied, keeping k x = fixed. The lowest four Landau bands are shown, while eight Landau bands were taken 
into account for the numerical calculation. The boxes depict the energy interval chosen for the wave packet dynamics in Fig. El 
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that leads to the same spectral and dynamical phenom- 
ena. To this end, we will study a limit where the classi- 
cal dynamics is not chaotic and where we will argue that 
avoided crossings are due to tunneling. 

IV. THE LIMIT OF NEGLIGIBLE LANDAU 
BAND COUPLING 



FIG. 2: Spreading of initially localized wave packets for 
magnetic modulation strengths (a) b x + b y — 1.1 and (b) 
b x + by = 7.0 at a fixed ratio b x /b y = 0.8 and v x = v y = 0. 
The wave packets are constructed from the eigenstates in the 
energy intervals depicted in Fig. U 



For small modulations strengths b x ,b y ,v x ,v y <C 1 we 
can neglect the coupling between Landau bands. Thus 
we ignore matrix elements with different v and v', leading 
to a discrete ID Schrodinger equation with next-nearest 
neighbor coupling for every Landau band, 



J 



tx{a n -x + a n+ i) + t2{a 

n—2 H~ d n -\-2) ~r 



21 i co> ( 27r-j^(n + k y )j + 2U 2 cos Un^-(n + k y ] 



a n — Ea TI 



(17) 



r 



with a n ^ 



a n - The parameters Ui, U%, t\, and strengths, 



t2 are given by the magnetic and electrostatic modulation 



XJ\ = b x si + v x s 2 
t\ = b y si + v y s 2 
U 2 - -b 2 x s 3 



(18) 
(19) 
(20) 
(21) 
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and E = E/hw c - [(1/ + |) + ($/$o) + ^) /(8tt)] • which solely depend on the Landau level index v and the 
Here we have introduced the following coefficients, magnetic flux ^/^o- 



-l 



(7T$ /$) l/ - 1 - fe 



+ 1 ^ — — In the special case t 2 = V% = 0, i.e. without magnetic 

k=Q (y — k)\{v + 1 — k)\ modulation, Eq. I|17H corresponds to the Harper equation 

" v\ (tt$ /$) l/ - k 13 • I n the general case, an important property of the 

s 2 = — — — — (23) Harper equation carries over: the Aubry duality I n 

k~o kl ^ ~ k > 1 order to see this, we multiply Eq. (UJ by e 2 ™(™+WV*o/* 

3 ^ - (- 47 r$ /$)-fc and sum over »• With 6 ™ = £n=-oc e 2 «(»+^)™^/* 

S3 = e aTr * > — — - , (24) this leads to a dual equation 

^ k\ [y — k)\ A 



U 2 (b n 



2ti cos I 2n — (m 



2f2 cos ( 47T— (m + fc x ) 



(25) 



where C/i , U 2 are exchanged with t±, t 2 . The selfdual 
point is given by U\ = t\, U2 = t%, i.e. when the mod- 
ulation strengths are equal in the x and y directions, 
Vx = Vy and b x = b y . Like in the Harper equation one 
expects here a fractal spectrum and fractal eigenfunc- 
tions, as visualized in Fig. The duality implies that 
from an eigenstate at one parameter set one can directly 
construct the eigenstate at the dual parameter set. In 
particular, a localized eigenstate leads to a dual state 
that is extended [34l |. 

An important property of the Harper equation, that for 
U\jt\ < 1 all eigenstates are extended and for U 2 /t 2 > 1 
all eigenstates are localized (at least for most irrational 
<&/$o) does not carry over. In fact, we see in Fig. 
that for a given parameter set localized and extended 
eigenstates coexist. This is accompanied by avoided band 
crossings and thus can be understood on the basis of the 
analysis in Ref. (l2T |. 

Two questions arise at this point: Firstly, one might 
wonder if in the special case of a purely magnetic modu- 
lation, like for a purely electrostatic modulation (Harper 
equation) , all eigenstates for a given parameter set are ei- 
ther localized or extended. The answer is given by Fig.Q] 
which shows an energy spectrum for the case of purely 



magnetic modulation, i.e. (U 2 /t 2 ) — {U\/ti) 2 . We still 
observe avoided band crossings and the coexistence of 
localized and extended states. We conclude that the ex- 
istence of the terms t 2 and U 2 in Eq. I|17J1 leads to this 
deviation from the qualitative properties of the Harper 
equation. As an aside, we note that we have to zoom into 
finer details of the spectrum to see this effect in the case 
of a purely magnetic modulation. 

The second and more important question stems from 
the fact that discrete ID Schrodinger equations like 
Eq. ifTTIl have a classical limit that is integrable (see be- 
low) , while it was previously argued (l2l H"^ that avoided 
band crossings are induced by a classically chaotic limit: 
Why do we observe avoided band crossings in this clas- 
sically non-chaotic case? 



V. CLASSICAL LIMIT, TUNNELING, AND 
AVOIDED BAND CROSSINGS 

A ID discrete Schrodinger equation can be written in a 
Hamiltonian form: Equation i|17H with x = 2TT^-{n + k y ) 
and ip(x) = a n can be written as, 



t 1 [^{x-27r^-)+^(x + 2TT^-)] +t 2 [ip{x-A-K^) + ip(x + A-K-^)] +2[Uicos(x) + U 2 cos(2x)]i;{x) = Eip(x). (26) 
Now we can utilize the translation operators exp (±2Tt-£-d x ) and exp (±4Tr^d x ) , yielding 



$0 3>o 
t\ cos(— i2ir — d x ) + t 2 cos(— iAir — d x ) + XJ\ cos(x) + U 2 cos(2x) 
<j> 4> 



ij)(x) = Eip(x). 



(27) 



r 



Comparing this with the time-independent Schrodinger with a momentum operator p — —ih c gd x and an effec- 

equation gives the Hamiltonian, tive Planck constant H e g — 2n^-. The corresponding 

H = 2[ti cos(p) + t 2 cos(2p) + Ui cos(x) + U 2 cos(2ir)] , 

(28) 
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Ui under variation of k x 
1.5, and Ui = 1.5 (a) and 
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n 

FIG. 3: Spectrum of Eq. fl3 vs 
with g/p = 377/610 for ti = 1, t 2 = 
magnification of the boxed region (b). At the selfdual point 
Ui = 1 all eigenstates are fractal. Many avoided band cross- 
ings can be seen and there are parameters values for which ex- 
tended (bottom left) and localized eigenstates (bottom right) 
coexist. 



classical Hamiltonian is one-dimensional and therefore 
integrable. The ubiquitous appearance of avoided band 
crossings, as in chaotic systems, is thus not expected. 

Still, we observe in Figs. and 21 some avoided band 
crossings. We explain this by studying the classical 
Hamiltonian Eq. (2HJ- For the case of a purely electric 
modulation (Harper equation), t% = U2 = 0, it has a sin- 
gle maximum and a single minimum within a unit cell. If, 
however, t% and U2 become non-zero, additional minima 
or maxima appear and classically energetic degeneracies 
become possible, as depicted in Fig. Quantum me- 
chanically these are split by tunneling and lead to avoided 
crossings under the variation of a parameter of the Hamil- 
tonian. The width of such avoided crossings, of course, 
depends on the effective Planck constant h c s- Here we 
have studied values of H c s, which are of the order of 1, ex- 
plaining the easily observable avoided crossings of Figs.EI 
and 01 From the classical Hamiltonian Eq. 11280 we see 
that a magnetic modulation, which leads to £2,^2 7^ 0, 
is a necessary condition for the appearance of classical 
energy degeneracies and quantum avoided crossings. We 




FIG. 4: Spectrum of Eq. I|17H vs. t 2 for a purely magnetic 
modulation, (U\/t\) 2 = under variation of k x with 

q/p = 377/610. Further parameters are t\ = 1 and U\ = 1.2. 
The inset shows avoided band crossings and the coexistence 
of extended and localized eigenstates. 




FIG. 5: Iso-potential lines of the classical Hamiltionian 
Eq. (EHJ with (a) fa = U 2 = and (b) t 2 ,U 2 + 0. While 
in (a) there are no energy degeneracies, in (b) several max- 
ima and minima occur and energy degeneracies exist, e.g. the 
one depicted by thick lines. 



note, that if one goes beyond the lowest Fourier compo- 
nents of the spatial modulation, Eqs. ltl5|l and ltl6ll . a 
purely electric modulation may also lead to energy de- 
generacies. 

Now the usual reasoning applies namely that 

avoided band crossings lead to metal-insulator transitions 
and the coexistence of localized and extended states, ex- 
plaining the observations in Figs. and 0] The dynamics 
of wave packets would show the corresponding changes 
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in transport direction similar to Fig. 



VI. CONCLUSION 

We derive a matrix equation for the energy spectrum 
of electrons in 2D periodic electric and magnetic fields. 
We observe avoided band crossings, metal-insulator tran- 
sitions, and changes in the spreading direction of wave 
packets. These phenomena, previously studied in the 
case without magnetic modulation, are now found in two 
new situations. 

Firstly, we find that they occur in a purely magnetic 
modulation. The origin being the classically chaotic dy- 
namics, which leads to avoided band crossings. Secondly, 
we study the limit of weak electric and magnetic mod- 
ulation, where the Landau band coupling can be ne- 
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